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We report on recent progress on perturbative QCD calculations of certain exclusive 
rare weak B meson decays involving hard photons. In the limit of a photon energy 
E-f much larger than Aqcd^ the amplitudes for such processes can be analyzed 
in a twist expansion in powers of A/E-y. The leading twist amplitude is given 
by the convolution of a hard scattering amplitude with the B meson light-cone 
wavefunction. This approach is applied to a calculation of the leptonic radiative 
B — ► •yii'i formfactors and to an estimate of the weak annihilation contribution to 
the penguin decays B — > p{u))^. As an application we discuss a few methods for 
constraining the unitarity triangle with exclusive radiative B decays. 

1 Introduction 

Rare radiative decays of B hadrons have been extensively studied as a possible 
way of probing possible new physics effects. Most effort has been concentrated 
on the theoretically cleaner inclusive decays, which can be analyzed in an f /mt 
expansion with the help of an operator product expansion. 

On the other hand, exclusive rare B decays are sensitive to long-distance 
QCD effects and the relevant amplitudes depend on details of the hadronic 
bound states. This makes their theoretical description considerably more 
model-dependent, and a full description is still lacking. 

In the following we describe a systematic treatment of exclusive decays in- 
volving one very energetic on-shell photon emitted from an internal quark line. 
In this situation the internal light quark moves very close to the light cone, 
and an expansion in powers of the small parameter A/Ej becomes usefuLl. 
This is analogous to the twist expansion for exclusive processes involving only 
light quarksB. 

The simplest process which can be studied along these lines is the leptonic 
radiative decayO B "yii'e, which is discussed in Sec. 2. A basic ingredient 
of the method is the light-cone description of a heavy-light meson, which is 
introduced in Sec. 2.1. In Sec. 2.2 we argue that the general structure of 
the leading-twist B —> jtiyi form factors is given by a convolution of the B 
light-cone wave function with a hard scattering amplitude Th- „ 

In Sec. 3 we describe an application of this method to the calculatiora 
of the weak annihilation contribution to the weak radiative decay B — > pj. 
Taking into account all possible long distance amplitudes compatible with 
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SU(3) symmetry, we present in Sec. 3.1 a few methods for constraitung the 
unitarity triangle using isospin violating effects in i? — s- p{iLj)j decayfu. 

2 Leading twist calculation of the B ^/Iv^ decay 
2.1 Heavy mesons light- cone wavefunctions 

The most general form for the Bethe-Salpeter wave function of a heavy-light 
meson B{v) can be written in terms of 4 scalar functions il^i{v-z, z^) {z = x—y) 

Mz,v) = {0\Tb4x)qp{y)\B{v)) = ((^i - yV)V2 + #3 + /V'4) 

^ J Q/3 

_ (1) 

In the heavy quark limit the structure of the wavefunction is further con- 
strained by the condition iZ-ip = ip. This reduces the number of independent 
structures to 2, which can be taken as V'1,2 

^Q/3(^,t')= {^(V-l + iZ-yCf •^)]V'2)75| ^. (2) 

The momentum space B meson wavefunction is defined as 

4,^p{k) = J d^ze''-'4^^p{v,z). (3) 

In the physical applications to be discussed in the following, this wavefunction 
is convoluted with a scattering amplitude T^f(fc) which does not depend on 
one light-cone component of the relative momentum in the bound state 
(e.g. Therefore, this component can be integrated over, which effectively 
puts z on the light cone (z+ — 0). We defined here light-cone coordinates as 
k± = k^ zLk^, which can be projected out by dotting into the light-cone basis 
vectors n = (1, 0, 0, 1) and h = (1, 0, 0, —1). This gives k+ = n ■ k,k- = n ■ k. 

It is convenient to write the wavefunction (||) as a linear combination of 
light-cone projectors in spinor space A+ = ^t;/^, A_ = j-r^i^ 

V'q/5(z+ -0,z_,zx -0) = |ii^(A+V^+(z_)+A_V7_(z_))75| ^. (4) 

with ^±{z-) — ipi ^ ^z_^2- [When used to compute a physical amplitude, 
only one of these functions will contribute to leading twist (e.g. This 
gives the most general expression for the B jMavefunction in the heavy quark 
limit, which is often found in the literature ua 

1 + y' f / 1 
il}ai3{v ■ z,z^) ^ —- ~ (-0+ - V^-))75 > . (5) 
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However, in the following we will use the form (j^ which has a simple in- 
terpretation in the parton model. With the kinematics adopted above, the 
wavefunction ip^{k^) gives the probability to find the light quark in the B 
meson carrying momentum Its moments are related to matrix elements 
of local operators 

{k'^)= / dk+k'^^j+{k+) = {0\qA+i/j5{in ■ DfK\B{v)) . (6) 
Jo 

The first two moments can be expressed in terms of known B hadronic 
parameters: {k^) = Jb^b and {k+) = ^AfBrUB, with A = — mb ~ 350 
MeV the binding energy of the b quark in a _B meson. The average of k^^ will 
play an important role in the following. Although it cannot be related to the 
matrix element of a local operator, it is possible to give a model-independent 
lower bound on its magnitude. Using the normalization conditions for N = 
0, 1 and the positivity of the distribution function one findsci 

For more specific (but less model-independent) predictions, some model has 
to be adopted for the distribution function We used an ansatz 

inspired by a quark model with harmonic oscillator potential = 
Nk+ exp (— 2Zr^(^+ ~ ^)'^)- '^^^'^ width parameter uj is varied as a; = 0.1 — 0.3 
MeV and a is constrained from the normalization conditions. Using A = 
0.3 - 0.4 MeV gives a = 0.05 - 0.5 MeV. 

2.2 Leading twist form factors in B ^ 'jii^i 

The simplest hard photon process which can be analyzed using an expansion 
in l/E^ is the radiative leptonic decay B ^iv^. This proceeds through the 
weak annihilation of the spectator quark in the B meson, as depicted in the 
quark diagrams in Fig. ^. ^ 

Interest in this decay was sparked by the observationEI that its branching 
ratio is enhanced relative to that for the leptonic decay B — > Hu^ with £ = e, ^. 
Although adding one photon introduces a factor of a = in the rate, it 
also removes the helicity suppression factor {m^/rriB)'^', the overall effect ia 
an enhancement in the rate of the leptonic radiative mode. Later workt3 
investigated this decay using a variety of approaches, both in connection with 
the prospect of constraining and/or \Vub\ and in the context of the weak 
radiative decays B p(w)7. 



(7) 
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(a) 



(b) 



Figure 1. Tree level contributions to the hard scattering amplitude for B — > -fiui decay. The 
cross denotes the weak current = q7;j(l — 75)b- The dominant contribution comes from 
the photon being emitted from the light quark line (a) while the diagram (b) is suppressed 
by the inverse heavy quark mass l/rrii,. 



The amplitude for B 
fv,A[E^) defined as 



1 



\/47ra 
1 



"fiv^ is parameterized by two formfactors 
{^{q,e)\q^^,h\B{v)) ^ie^,^pse*''v%' fv{E^) (8) 



\/47ra 



{l{q.e)\q-t^-t^h\B{v)) = [q^{ve*)-e*^{vq)\fA{E^). 



(9) 



Computing the contribution of the two diagrams in Fig. I with the B wave- 
function one finds for the formfactors (^), (^) 



fv{E,)^fA{E,)^l^\Q,R- 



2Ey \ * TOf, 



(10) 



where Qq,Qb are the spectator- and heavy-quark electric charges and the 
hadronic parameter R describes the photon coupling to the light quark. This 
is given by a convolution of the leading twist B meson wavefunction with a 
hard-scattering amplitude Th{E^, k-^-) 



dk+ — ; TH{Ey,k+) 

k+ 



(11) 



The hard scattering amplitude Th{Ej, fc+) is given to one- loop order b; 



TH{E^,k+) = 1 



4Tr 



2E^ 



/2K 



-log^(^)+ilog(^if )-^) • (12) 



47r" 
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Using the model distribution function described in Sec. 2.1 one finds at tree 
level R{E^) = 2-4 GeV-\ while the lower bound (0) predicts i? > 2.1 
GeV-i (corresponding to A = 350 MeV). 

We briefly discuss in the following a few interesting consequences of these 
results. 

a) To leading twist the fornifactors fv,A{Ej) scale like with a pro- 
portionality coefficient which is independent on the heavy quark flavor (up to 
l/rrib corrections). This property can be used to determine the CKM matrix 
element. \ Vub\ from a comparison of the B —f jivi and D — > jlvi photon 
spectrao. 

b) The form factors of the vector and axial current are equal to leading 
twist. This is in contrast to their behaviour in the low Ej region, where 
they receive contributions from intermediate states with different quantum 
numbersO (J^ = 1~ for fv and = 1+ for /a)- 

The equality of the form factors (|l^ is a particular case of a symmetry 
relation analomus to those discussed for the soft components of semileptonic 
formfactors inllll. The gluon couplings of a quark moving close to the light 
cone possess a higher symmctr5c3. This can be formalized by going over to 
an effeclive theory, which should include in addition to the soft gluon modes 
(LEETIIj), also collinear gluons. A_-|Complete discussion including collinear 
gluons has been given only recentlyO. 

There is however an important difference between the status of the sym- 
metry relation fy = /a among the leptonic radiative formfactors^-and the 
analogous symmetry relations for the soft semileptonic formfactorsllll. While 
the latter receive symmetry breaking corrections from hard gluon exchange 
(which have been computed recentlycl), the former relation is not changed by 
such effects, as checked by explicit calculation to one-loop ordcrQ. 

Hadronic amplitudes similar to R in ( pT| ) appear in many physical quanti- 
ties involving long-distance effects iiiduced by hard photon or gluon emission 
from weak annihilation topologiesEil. [When computed in the quark model, 
this quantity appears as the inverse of the constituent quark mass R — > l/m^.] 
In the following section we study such an important application, to long- 
distance effects in exclusive penguin induced decays b d^y. 

3 Long-distance contributions to the B pj decay 

An important class of weak radiative decays are those mediated by the penguin 
mechanism b — > s{d)^ 

= =^VtbV^:Cr^F^'•'sa^,PRb + Hid. , (13) 
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q 



(d) (e) 

Figure 2. Long distance contributions to B ^ V-y weak radiative decays. The crosses 
denote the possible attachments of the photon to the internal quark lines. The quark 
diagrams are denoted in text as a) weak annihilation (A); b) W-exchange (E); c) penguin 
amplitudes containing u-quark (Pu^) and c-quark (Pi*^) loops. We distinguish between 
amplitudes with the photon attaching to the spectator quark (j = 1) and to the quark in 
the loop (i = 2). d) annihilation penguin amplitudes {PAu,c); e) gluon penguin amplitudes 
{M(*)). 



with small long-distance contributions expected from time-ordered products of 
the weak nonlcptonic Hamiltonian with the minimal electromagnetic coupling 
(denoted as Tii.d.)- 

We show in Fig. ^ the different quark diagrams responsible for long- 
distance contributions to a typical weak radiative decay B Vj. InJhe 
SU(3) limit any such amplitude can be written as a linear combinationO of 
the amplitudes in Fig. ^ with CKM factors Vi (for each photon helicity A) 

A(B^F7a)= E ^i'^K E C,,gM,,gX. (14) 

q=u,c,t Mi.q=A,E,P,,PA, 

This is analogous to a similar jdficomposition of nonlcptonic B decay ampli- 
tudes into graphic amplitudeaij, and is equivalent to a more usual SU(3) 
analysis in terms of reduced matrix elements. 

The long-distance amplitudes appearing in ( p^ are notoriously diffi- 
cult to calculate. JX^s^i have been estimated JJSing various methods_,such 
as QCD sum rplcalj, vector meson dominancetll, perturbative QCDE3 and 
Regge methodsE^I. A few sample results are tabulated in Table 1, separately 
for the two helicities of the photon X = L,R. These calculations show that 
the dominant long-distance amplitude in 6 — s- o?7 decays comes from the weak 



proc: submitted to World Scientific on February 7, 2008 6 



Table 1. Estimates of the short-distance and long-distance amplitudes in B ~* p-y decays 
(in units of 10~^ GeV). The estimates of the WA and VF-exchange amplitudes Ax and E\ 
used R = 2.5 GeV~^. The penguin type ampUtudes Pu and Pc have been estimated using 
vector meson dominance. 



Photon hclicity 






\Pux\ 




\Ex\ 


X = L 


1.8 


0.16 


0.03 


0.6 


0.05 







0.04 


0.007 


0.07 


0.007 



annihilation graph A in Fig. ]^{a). 

It is therefore rather fortunate that the weak annihilation amplitude can 
be computed in an esenlially model-independent way. In the factorization 
approximation, one findsO 



G C 



E^{fv±fA)) 



(15) 



and 



with fv.A{E^) the radiative leptonic form factors defined in (|8|), 
C'i(mb) = —0.29, Ci(mb) = 1.13 are Wilson coefficients in the weak nonlep- 
tonic Hamiltonian. Nonfactorizable corrections to this result arise from hard 
gluoiis connecting the initial and final state quarks, but appear only at higher 
twisifi. Furthermore, to leading order of a twist expansion for the radiative 
leptonic formfactors fv.A (|lO|), the weak annihilation amplitude couples pre- 
dominantly to left-handed photons. A similar suppression of the right-handed 
helicity amplitudes is noted for all other long-distance contributions (see Ta- 
ble 1), which has implications for propop.ls to search for new physics through 
photon hclicity effects in 6 ^ 57 decayaHj. 

The main motivation for measuring exclusive weak B radiative decays is 
connected with the possibility of extracting the CKM matrix element |Vtt;|. 
The weak annihilation amplitude A introduces the most significant theoretical 
uncertainty in such a determination. Keeping only the leading long-distance 
amplitude, one finds for the ratio of charge-averaged rates 



td 



Vt, 



R 



SU{3) 



V,hV 



b^ud 



VtdVg 



cos a cos 



cl>A + 0{e\) 



(16) 

with EAe*"^-^ = Al/Pil and Rsup^ = 0.76±0.22 an SU(3) breaking parameter 
in the short-distance amplitudecil. Using the estimates of Table 1 one finds 
ea — 0.3, which gives an uncertainty of about 15% in this determination of 
\Vtd\ from long-distance effects. The amplitude A induces also a direct CP 
asymmetry in charged B decay Ac p — 2 X/''Y/t'^ £a sin a sin (/)a , which can 
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be as large as 30% for optimal values of the weak and strong phases a, (f)A- 



3. 1 Constraining the CKM matrix with exclusive weak radiative B decays 

Interference between the e.m. penguin and long-distance amplitudes can pro- 
duce isospin breaking effects in ~^ p7 decays. Some of the latter contribute 
with a different weak phase than the former, which led to the suggestion to 
use isospin breaking in these decays in order to extract information about 
CKM parametersB. 

A more careful treatmentcl shows that such an approach could receive con- 
taminations from additional long-distance effects with the same weak phase 
as the short-distance amplitude. These amplitudes can be related by SU(3) 
symmetry to isospin-breaking effects in i? ^ K*^ decays. Eimerimentaldata 
on these modes have been recently reported by the CLE0c3, BaBaiO and 
Bella^j collaborations, which find 

B{B^ K*^j) = (3.761°:^^ ± 0.28) x 10"^ [CLEO) (17) 

(2.87 ± LiOtoTo) X 10"^ {BELLE) 

B{B° K*"-f) = (4.55^;°:^^ ± 0.34) x lO^^ (CLEO) (18) 

(4.94 ± 0.93toil) X 10"^ (BELLE) 

(5.2 ± 0.82 ± 0.47) x 10^^ (BABAR) 

The long-distance amplitudes responsible for the difference in rate between 
charged and neutral B —>■ K*^ modes are shown in Figs, ^(c), (e), and contain 
a charm loop or a gluon penguin in which the photon attaches to the spectator 
quark. The preliminary data (p^, ( p^ indicate that these effects could be 
significant. It has been proposecB therefore to eliminate them by combining 
B fry with B K*^ data, by forming the combined ratio 

^ B{B^ ^ p^7) B{B^ ^ K*°7) 

~ 2B{B^ p^-i) B{B± if *±7) ' 



1 - 2epA^ — 2e^ 



VtbV, 



td 



cos a cos (f)A + 0{e1) . 



We expanded here to linear order in the ratios of long-/short-distance am- 
plitudes Si. The residual contamination from the (OZI-suppressed) penguin- 
annihilation amplitude £pA^ Fig. ^(d) can be expected to be very small. An 
upper bound on its size can be given in terms of experimental data on Bg 
decays as jepyij^ < 2r(B^ ^ p°-i)/T(B^ ~* K*^-f). 

Any measurement of the ratio ( |l9| ) different from 1 can be translated into 
a constraint on the CKM factors in the last term. The cleanest approach 
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0.2 0.4 0.6 0.8 1.0 p 

Figure 3. Illustrative example of a constraint on the CKM parameters from isospin breaking 
in B — > p{u))^ decays. The shaded area is the 68% CL contour obtained from a global fit 
of the unitarity triangle^^. The region contained between the dotted lines is the exclusion 
region corresponding to \Rp — 1\/£a = 0.1. 



involves extracting the factor 
^ivi and B ^ K*j decays 



ea from a combination of data in _B 

This resuhs into a bound on the weak phase a 
which excludes values around a = 90°. 

A simpler (although less model-independent) approach would combine 
theoretical calculations of ea using , in order to constrain the combination 

of CKM parameters 



cos a. A typical region in the (p, rf) plane which 

can be excluded in this way is shown in Fig. ^ This largely overlaps with 
the area presently favored by global fits of the unitarity triangle, and can 
therefore be expected to be useful in further constraining it, as more data on 
weak radiative B decays become available. 
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